Confinement vs Deconfinement of Cooper Pairs in One-Dimensional Spin-3/2 

Fermionic Cold Atoms 
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The phase diagram of spin-3/2 fermionic cold atoms trapped in a one-dimensional optical lattice 
is investigated at quarter filling (one atom per site) by means of large-scale numerical simulations. 
In full agreement with a recent low-energy approach, we find two phases with confined and decon- 
fined Cooper pairs separated by an Ising quantum phase transition. The leading instability in the 
confined phase is an atomic-density wave with subdominant quartet superfluid instability made of 
four fermions. Finally, we reveal the existence of a bond-ordered Mott insulating phase in some part 
of the repulsive regime. 
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Loading cold atomic gases into optical lattice allows 
for the realization of bosonic and fermionic lattice models 
and the experimental study of exotic quantum phases 1 -. 
Ultracold atomic systems also offer an opportunity to in- 
vestigate the effect of spin degeneracy since the atomic 
total angular momentum F can be larger than 1/2 re- 
sulting in 2F+ 1 hyperfine states. This high-spin physics 
is expected to stabilize novel exotic phases. In this re- 
spect, various superfluid condensates, Mott insulating 
phases, and interesting vortex structures have been found 
in spinor bosonic atoms with F > l 2,3 . These theoreti- 
cal predictions might be checked in the context of Bose- 
Einstein condensates of sodium, rubidium atoms and in 
spin-3 atom of 52 Cr— . The spin-degeneracy in fermionic 
atoms is also expected to give rise to some interesting 
superfluid and Mott phases. In particular, a molecular 
superfluid phase might be stabilized where more than 
two fermions form a bound state. Though such non- 
trivial superfluid behavior has been previously found in 
different contexts^, it has been advocated recently that 
the formation of bound-state of Cooper pairs is likely to 
occur in general half-integer F > 1/2 ultracold atomic 
fermionic systemsSi^. In the spin F = 3/2 case, it has 
been predicted on the basis of a low-energy study&£ in 
one dimension that a quartetting superfluid phase, i.e. 
a bound-state of two Cooper pairs, might be stabilized 
by strong enough attractive interactions. The simplest 
lattice Hamiltonian to describe spin-3/2 atoms with s- 
wave scattering interactions in ID optical lattice takes 
the form of a Hubbard- like model 2 : 
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where c\ i is the fermion creation operator corresponding 



to the four hyperfine states a = ±1/2, ±3/2. The sin- 
glet and quintet pairing operators in Eq. ([TJ) are defined 
through the Clebsch-Gordan coefficient for two indistin- 
guishable particles: P] M% = J2 a p{JM\F, F; a/^c^cj^. 
As it appears, it is more enlightening to express model ([1]) 
in terms of the density (jii = J2 a c a i°a.i) an d the singlet 
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with U = 2U 2 and V = U - U 2 . Model © generi- 
cally displays an exact SO(5) extended spin symmetry 
and an SU(4) symmetry in the particular case Uq = Ui, 
i.e. V = 0^. In sharp contrast with the spin F = 1/2 
case where both interacting terms in Eq. © are pro- 
portional, these terms are independent for F = 3/2 
and strongly compete. While the V^-term favors the 
pairing of two fermions for negative V, an attractive 
[/-interaction might favor the formation of a quartet 
Qi = c_3 i c_i 4 Ci ^03 4 . In fact, it has been recognized in 
Ref. l6| that the above competition reveals itself through 
a non-trivial discrete symmetry of the problem. Indeed, 
model ([2]) possesses, on top of the SO(5) symmetry, a 
Z2 discrete symmetry U: c a .i — > e l7r / 2 c Q .i which plays a 
crucial role in the low-energy physics since, as Pj is odd 
under U, the formation of a quasi-long range BCS phase 
requires hi to be spontaneously broken. When U is un- 
broken the BCS instability is strongly suppressed and the 
leading superfluid instability is made of four fermions i.e. 
a quartet which is even under hi. Such a two-phase struc- 
ture has been recently predicted away from half-filling in 
the weak coupling limit by means of a low-energy ap- 
proac h 6 i 7 i 10 . 
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FIG. 1: (Color online) Phase diagram of the spin-3/2 Hub- 
bard chain ([2]) at quarter filling from QMC and DMRG cal- 
culations for V < (see text for definitions). 
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FIG. 2: (Color online) Correlation functions for U = — 1.5t 
and V — obtained by DMRG. (a) Power-law behavior of the 
quartet correlations parametrized by Luttinger exponent Kd- 
(b) Short-range behavior of the one-particle Green function 
and pairing correlations (£ denotes the correlation lengths). 



In this letter, we investigate numerically the phase 
diagram of model ^ for V < at quarter filling 
(one atom per site) by means of quantum Monte-Carlo 
(QMC) and Density-Matrix Renormalization Group 11 
(DMRG) simulations. Physical properties are investi- 
gated by computing the one-particle, density, pairing 
as well as quartetting correlation functions, respectively: 
G{x) = (cl lCa , l+x ), N(x) = ( ni n l+x ), P(x) = {PiPj +x ) 

and Q(x) = (QiQ\ +x )- For the QMC simulations, we 
used the projector auxiliary field QMC algorithm (see 
Ref. Oj] for the details of the algorithm) in the regime 
V < and U < — 3 V/ 4 where the fermionic algo- 
rithm has no sign problem 9 . We have studied periodic 
chains with linear size up to L = 180 with a typical 
projection parameter Qt = 10 and a Trotter time slice 
At = 0.05. Most of DMRG calculations were performed 
on open chains with L = 60 sites and keeping M = 1400 
states^. The resulting phase diagram at quarter filling 
is presented in Fig. [T] and we now turn to the discussion 
of the physical properties of the different phases. 

Confined Phase - The phase with U < and small 
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FIG. 3: (a) Fourier transform N(q) of the density correla- 
tions obtained from QMC (U = —t and V = 0). The linear 
dispersion at small q gives access to the Luttinger parameter 
Kd- Insert : the scaling of the peak at 2&f vs L signals an 
ADW phase, (b) Luttinger exponent Kd as a function as U. 
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FIG. 4: (Color online) Correlation functions for [7 = and 
V = — 3t from DMRG. (a) Pair and quartet correlations are 
algebraic with Q(x) ~ P(x) 4 . (b) Short-range behavior of 
the one-particle Green function. 



\V I is characterized by the existence of a spin gap and an 
unbroken Z2 discrete symmetry U which marks the onset 
of an Atomic Density Wave (ADW) and quartet super- 
fluid quasi-long-range orderings. Indeed, for the typical 
value of U = —1.5t and V = 0, we observe in Fig. [5{b) 
that both the pairing correlations and the one-particle 
Green function decay exponentially with distance. In 
contrast, the quartet correlations are algebraic as it can 
be seen in Fig. [DJa). We have checked by a direct evalua- 
tion that the four-particle gap vanishes. We can thus 
deduce that the short-range character of P(x) is due 
to the confinement of Cooper pairs which stems from 
the unbroken U symmetry. The above results extend in 
the whole confined phase (squares in Fig. [T]). In this 
phase, the superfluid instability is of a molecular type 
made of four fermions: a quartet. However, the den- 
sity correlations also display a power-law behavior with 
dominant oscillations at 2k p — "ft/2 as it is clearly seen 
from the Fourier transform N(q) of N(x) presented in 
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FIG. 5: (a) Fourier transform of the density correlations ob- 
tained by QMC for U = and V = 1.5*. (b) The Luttinger 
parameter Kd as a function of V when U — 0. 



Fig. [3ja). The question that naturally arises is which in- 
stability dominates in this phase. The answer depends 
on the value of the non-universal Luttinger parameter 
Kd which stems from the critical behavior of the den- 
sity degrees of freedom. Indeed, the quartet and 2kp- 
ADW equal-time correlation functions have been found 
in Refs£^ to behave at long distance as Q(x) ~ x~ 2 ' Kd 
and Nix) ~ cos(Trx/2)x~ Kd ^ 2 . Therefore a quartet su- 
perfluid phase with dominant quartet correlations re- 
quires Kd > 2. The value of Kd has been computed 
in QMC using the formula: 



i^lim^M. 

4 9-o q ■ 



(3) 



where the factor 4 comes from the four spin states. This 
procedure has been shown to be very accurate for the 
spin- 1/2 Hubbard modeii^. For DMRG calculations, Kd 
can be independently obtained from the power-law be- 
havior of the quartet correlations Q(x). Both QMC and 
DMRG results are shown for example on the SU(4) in- 
variant line (V = 0) in Fig. |3fb). We find that QMC 
works better than DMRG in the weak coupling limit and 
is in excellent agreement with the perturbative estimate: 
K^ 2 =l + [V + 3C/]/(\/27rf)£il. For larger \U\, DMRG 
is more accurate and we found that Kd saturates at the 
value Kd — 1.6. Note that the perturbative estimate 
fails beyond \U\ > t so that numerical approaches be- 
come necessary to estimate Kd- For V ^ 0, Kd also 
saturates at strong couplings to values smaller than 2. 
We therefore conclude that, though the quartet correla- 
tions are quasi-long ranged, the dominant instability in 
the confined phase is a 2fcp-ADW. 

Deconfined Phase - By allowing V to be sufficiently 
negative, one can enter a second phase where the 
one-particle gap is still finite (see Fig. 0Jb) for U = 
0, V = —3t) but the two-particle gap vanishes. In this 
phase, pairing correlations become algebraic with P{x) ~ 
x -x/2K d ag s j 10wn on pig gja) and Q(x) remains critical 
with Q(x) ~ P(x) 4 which is the prediction of the low- 
energy approach. In contrast to the ADW phase, there is 
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FIG. 6: (Color online) The BCS-ADW transition from DMRG 
computations along the V — —2t line, (a) Normalized ratio 
R(x) — P(x) 4 /Q(x) displaying a critical behavior at the tran- 
sition, (b) In the bulk (at site x = 45), R(x) is proportional 
to U - U c for U C <U <U* with U* = -t and U c = -1.191 
(c) Plot of R(U) = a{U) s vs t = ((7- U C )/(U* - U c ) where 
a is the Ising order parameter. 



no diverging signal at 2kp = tt/2 in N(q) (see Fig.(5ja)). 
We thus conclude that there is still a spin-gap and the Z2 
symmetry hi is now spontaneously broken which leads to 
the formation of a quasi-long-range BCS phase. In addi- 
tion, there is also an ADW instability at 4kp = w (see 
Fig. EJa) where N(q) has a maximum at q = 7r) which 
has a power-law decay N(x) ~ cos(7ra;)x~ 2 - ft - ti . We thus 
need to compute numerically Kd to fully characterize the 
dominant instability of this phase. As in the previous 
phase, the Luttinger parameter Kd can be extracted ei- 
ther from Eq. ([3]) (QMC) or from pairing correlations 
(DMRG). As shown on Fig. Efb), both results are com- 
patible and agree with the perturbative estimate when 
|V| < t. We find that Kd > 1/2 so that the dominant 
instability in this phase is the BCS singlet pairing. 

Quantum phase transition - The striking feature of 
the phase diagram for attractive U, V interactions is the 
change of status of the Z2 symmetry IA which is sponta- 
neously broken (resp. unbroken) in the deconfined (resp. 
confined) phase. We thus expect a quantum phase tran- 
sition in the Ising universality class between these two 
phases. In fact it has been shown in Ref. that the 
order parameter o~(x) associated with the U symmetry, 
though being non-local in terms of the lattice fermions, 
can be extracted from the long-distance behavior of the 
ratio R(x) = P(x) A / Q(x). In the confined phase where U 
is unbroken, (<r[x)) — and R(x) ~ (a(x)a(O)} 4 decays 
exponentially with distance. In the deconfined phase, 
(ct(x)} = it / and R(x) ~ a 8 . Finally, it has been 
found in Ref. that at the transition the ratio displays 
an universal power-law behavior: R{x) ~ 1/x. We have 
computed numerically this ratio by DMRG for various 
parameters to determine the transition line in Fig. [T] 
The results of Fig. [H^a) clearly show an excellent agree- 
ment with the predictions of the low-energy approach. 
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FIG. 7: (Color online) Mott phase for V < 0. (a) By increas- 
ing U, pairing correlations obtained from DMRG change from 
algebraic to exponential decaying, (b-d) The BOW Mott 
transition is seen from DMRG computations with the appear- 
ance of the 4/cf order parameter t(x) across the transition. 



In particular, we observe that R{x) ~ 1/x near the crit- 
ical point. In the deconfined phase, R{x) saturates at 
large distance as it should and is almost independent of 
x (R(x) r~j a(U) s ) when one enters the critical regime 
(for U c < U < U* in Fig. [6(b)). The plot in Fig. E^c) 
demonstrates that <j{U) ~ (U — Uc) 1 / 8 in full agreement 
with Ising criticality. In this respect, the situation is in 
sharp contrast with the F — 1/2 well-known case where 
the 2fcf-ADW and BCS instabilities coexist for attractive 
interaction^. 

Mott phase - At quarter-filling, a Mott transition 
might take place if Kd < 1/2 with the formation of a 
density gap££ For the repulsive SU(4) Hubbard chain 



(V = 0), the QMC study of Ref. QJ found a transition 
from a gapless spin-density wave (SDW) to a generalized 
Mott SDW with three gapless spin modes (see Fig. [I]). 
For V < 0, we expect an entirely different Mott phase 
due to the presence of a spin gap and the breaking of 
the Z2 symmetry U. In the Mott region in Fig. [T] (full 
circles), the BCS singlet pairing becomes short-ranged 
as shown in Fig. [7(a) and we find that, as the density 
gap opens, the local density almost does not fluctuate 
and N(x) ~ 1. In contrast, the local kinetic bond, 

t( x ) = (J2 a c a ,x+i c o,x + H.c), orders with a Ak F = it 
modulation reminiscent of a doubly degenerate ground 
state as it can be seen in Fig. [7(b-d) for V = —2t. We 
therefore conclude on the emergence of a bond-ordering 
Mott phase with periodicity two (Mott BOW in Fig. [TJ. 

Concluding remarks - We conclude this letter in em- 
phasizing that the existence of a quartet superfluid phase 
where quartet correlations dominate over the 2kp-ADW 
instability relies on the non-universal Luttinger param- 
eter Kd- Though Kd < 2 at quarter filling, we expect 
that at sufficiently low densities, Kd may become larger 
than 2, which marks the onset of the quartetting phase. 
The formation of this exotic phase will be discussed in a 
forthcoming paper. 
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